Astrophysics > Earth and Planetary Astrophysics
[Submitted on 22 Aug 2022 (v1), last revised 30 Aug 2022 (this version, v2)]
Title:Trapping (sub-)Neptunes similar to TOI-216b at the inner disk rim: Implications for the disk viscosity and the Neptunian desert
View PDFAbstract:[Abridged] The occurrence rate of observed sub-Neptunes has a break at 0.1 au, which is often attributed to a migration trap at the inner rim of protoplanetary disks where a positive co-rotation torque prevents inward migration. We argue that conditions in inner disk regions are such that sub-Neptunes are likely to open gaps, lose the support of the co-rotation torque as their co-rotation regions become depleted, and the trapping efficiency then becomes uncertain. We study what it takes to trap such gap-opening planets at the inner disk rim. We performed 2D locally isothermal and non-isothermal hydrodynamic simulations of planet migration. A viscosity transition was introduced in the disk to (i) create a density drop and (ii) mimic the viscosity increase as the planet migrated from a dead zone towards a region with active magneto-rotational instability (MRI). We chose TOI-216b as a Neptune-like upper-limit test case, but we also explored different planetary masses, both on fixed and evolving orbits. For planet-to-star mass ratios $q\simeq(4$-$8)\times10^{-5}$, the density drop at the disk rim becomes reshaped due to a gap opening and is often replaced with a small density bump centered on the planet's corotation. Trapping is possible only if the bump retains enough gas mass and if the co-rotation region becomes azimuthally asymmetric, with an island of librating streamlines that accumulate a gas overdensity ahead of the planet. The overdensity exerts a positive torque that can counteract the negative torque of spiral arms. In our model, efficient trapping depends on the $\alpha$ viscosity and its contrast across the viscosity transition. In order to trap TOI-216b, $\alpha_{\mathrm{DZ}}=10^{-3}$ in the dead zone requires $\alpha_{\mathrm{MRI}}\gtrsim5\times10^{-2}$ in the MRI-active zone. If $\alpha_{\mathrm{DZ}}=5\times10^{-4}$, $\alpha_{\mathrm{MRI}}\gtrsim7.5\times10^{-2}$ is needed.
Submission history
From: Ondřej Chrenko [view email][v1] Mon, 22 Aug 2022 12:35:13 UTC (15,623 KB)
[v2] Tue, 30 Aug 2022 10:12:07 UTC (15,624 KB)
Current browse context:
astro-ph.EP
Change to browse by:
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
IArxiv Recommender
(What is IArxiv?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.